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Progress Report for NAGW-2306 


The Embedded Young Stars in the Taurus- Auriga Molecular Cloud: 

Models for Spectral Energy Distributions 

In this paper, which appeared in Kenyon et al. (1993a), we extended the pioneering 
work of Adams et al. (1987) to model the spectral energy distributions (SEDs) of all of the 
known Class I (protostellar) sources in the Taurus molecular cloud. The key ingredient of 
the models which produces reasonable fits rotation, which causes material to fall directly 
onto the disk, not the star, reducing the extinction in the inner envelope and allowing 
more short-wavelength light to escape. We adopted the Terebey, Shu, & Cassen (1984 = 
TSC) density distribution for a rotating, infalling envelope. At large distances, the TSC 
model is in essentially spherical steady free-fall, with a density distribution p za r~ 3 / 2 . On 
small scales rotation becomes important, and most of the material lands on a disk rather 
than directly onto the star. The centrifugal radius R c denotes the maximum radius 1o 
which infalling material lands on the disk at a given instant. The basic TSC model is 
defined by four parameters: the density at a reference level pi, the centrifugal radius R c , 
the total system luminosity L, and the inclination i of the polar (rotational) axis to the 
line of sight. 

We calculated the radiative equilibrium temperature distribution from the spherical 
average of the TSC density distribution (Adams & Shu 1986). The resulting spherically- 
symmetric temperature distribution then provides the source function to obtain the 
emergent spectrum at a given inclination angle i from the formal solution of the transfer 
equation, using the exact density (opacity) distribution, which is axisymmetric but not 
spherically-symmetric. We assume that the luminosity originates entirely in the central 
source, which radiates isotropically. To minimize the number of parameters we adopt 
the standard Draine & Lee (1984) dust opacities and assume that dust is destroyed at 
temperatures ^ 1600 K. 

Our results showed that the SEDs of the protostar candidates in Taurus can be 
reproduced with TSC models having infa.ll rates close to the values predicted by the 
theory of isothermal cloud collapse (e.g., Shu 1977). We also found that the cloud 
angular momenta had to be such as to produce infall to disks with R c ~ 100 A. II. 
This result is plausible, since the typical estimated sizes of T Tauri disks are of this 
order. However, the results indicated that most of the protostar candidates had large 
R c envelopes, which is not predicted by the TSC model. In addition, the model results 
formally indicated a tendency to observe Class I sources more pole-on, but this was 
probably due to our neglect of bipolar outflow cavities in the model, as suggested by 
scattered light calculations (Kenyon et al. 1993b). 
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Flat Spectrum T Tauri Stars: The Case for Infall 


We showed that the mid- to far-infrared fluxes of “flat spectrum” T Tauri stars 
can be explained by radiative equilibrium emission from infalling dusty envelopes. Infall 
eliminates the need for accretion disks with non-standard temperature distributions. The 
simplicity and power of this explanation indicates that models employing “active” disks, 
in which the temperature distribution is a parameterized power law, should be invoked 
with caution. Infall also naturally explains the scattered light nebulae detected around 
many flat-spectrum sources. To match the observed spectra, material must fall onto a 
disk rather than the central star, as expected for collapse of a rotating molecular cloud. 
It may be necessary to invoke cavities in the envelopes to explain the strength of optical 
and near-infrared emission; these cavities could be produced by the powerful bipolar 
outflows commonly observed from young stars. If viewed along the cavity, a source may 
be lightly extincted at visual wavelengths, while still accreting substantial amounts of 
material from the envelope. Infall may also be needed to explain the infrared-bright 
companions of many optical T Tauri stars. This picture suggests that many of the flat 
spectrum sources are “protostars” - young stellar objects surrounded by dusty infalling 
envelopes of substantial mass. 

Our results showed that the SEDs of the protostar candidates in Taurus can be 
reproduced with TSC models having infall rates close to the values predicted by the 
theory of isothermal cloud collapse (e.g., Shu 1977). We also found that the cloud 
angular momenta had to be such as to produce infall to disks with R c ~ 100 A.IJ. 
This result is plausible, since the typical estimated sizes of T Tauri disks are of this 
order. However, the results indicated that most of the protostar candidates had large 
R c envelopes, which is not predicted by the TSC model. In addition, the model results 
formally indicated a tendency to observe Class I sources more pole-on, but this was 
probably due to our neglect of bipolar outflow cavities in the model, as suggested by 
scattered light calculations (Kenyon et al. 1993b). 

This paper has been submitted to The Astrophysical Journal. 


Flattened Infalling Envelope Models 

Evidence is mounting that many young stars have dusty flattened envelopes around 
them, with sizes of a few thousand A.U. For example, the T Tauri star HL Tau has been 
found to have such an envelope in 13 CO radio interferometry. Sargent &; Beckwith (1991) 
interpreted this structure as a rotating flattened disk, but Hayashi et al. (1993) argued 
that the motions of this material are more consistent with infall than rotation. Hayashi 
et al. suggest that this structure is a “psuedodisk” of the type predicted by Galli & Shu 
( 1993a, b; see also Fiedler Sz Mouschovias 1993). In the Galli-Shu model, even if magnetic 
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fields are not of sufficient strength to prevent collapse, they can deflect infalling material 
into a “pseudodisk” which is not rotationally supported and then must continue falling 
in. 

Some further thought suggests that pseudodisks might be a common occurence in 
star formation even without the presence of magnetic fields. If collapse is initiated in 
a flattened, rather than spherical cloud, it is likely that the flattening will be enhanced 
during collapse rather than suppressed. That such flattening might be general is indicated 
by the studies of star-forming molecular cloud cores which by Myers et al. (1991), which 
indicated average aspect ratios ~ 2 : 1. Beyond this, the idea that stars form from clouds 
fragmenting from molecular cloud sheets or filaments (Larson 1985) has gained increasing 
interest and support (e.g., Bonnell & Bastien 1992). It is likely that a fragmenting sheet 
could also collapse to form a pseudodisk-type of infalling envelope. 

To examine this possibility further, we are collaborating with Dr. Alan Boss of 
the Carnegie Institution (DTM). Dr. Boss is calculating the collapse of a fragmenting 
sheet configuration with his time-dependent hydrodynamics code. For initial explorations 
isothermal collapse should be an adequate assuption. The dynamical models have been 
calculated to many free-fall times already for one set of initial conditions. An analytic 
approximation to the infall has been developed which reproduces many of the main 
features of the collapse models, and we are using this approximation to modify the 
radiative transfer solutions to investigate the qualitative effects expected. Ultimately we 
will use the density structures calculated by Dr. Boss at specific snapshots in time as 
the input for the radiative transfer models. The resulting predictions will be compared 
with observations in the optical, near-infrared, and far-infrared spectral regions. 

The collapse of the flat layer produces somewhat different results from the model of 
Terebey et al. (1984), which begins with an initially spherical cloud. We tend not to 
get an “inside-out” collapse initially, even though we have specified the location of the 
central mass by the assumption of axial symmetry. The outer layers react to the growth 
of the central perturbation before an expansion wave can arrive there because the gravity 
changes are felt immediately, unlike the spherical case. More importantly, the infall rate 
now depends not only on the sound speed (or gas temperature) but on the surface density 
as well, since for a given sound speed one may have a wide range of surface densities. 

This picture of a collapsing “thick disk” or toroid may also explain the infalling 
envelope of HL Tau on size scales ~ 1000 — 2000 A.U. The models also suggest that the 
infalling envelope will become increasingly flattened as the collapse proceeds (see Figure 
1). At long times the density distributions are reminiscent of the “flared disks” invoked 
by Kenyon Sz Hartmann (1987) to explain the far-infrared excesses of many T Tauri stars. 
It is conceivable that such relatively thick infalling structures could explain the infrared 
emission in excess of that expected from a flat, steady disk. 
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Figure 1 - density distribution in the meridional plane for the collapse model described 
in the text, at a time approximately 7 free-fall times after the start of the calculation. 
Each contour is a factor of two different than its neighbor, with the highest densities 
along the horizontal (equatorial) axis. The calculation shows the tendency for a flattened 
infalling envelope to develop, along with an axial hole in the envelope. 
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